Schròdinger Equation 
with the Potential V(r) = ar 2 + òr~ 4 + cr~ 6 
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Abstract 

By making use of an ansatz for the eigenfunction, we obtain the exact so- 
lutions to the Schròdinger equation with the anharmonic potential, V(r) = 
ar 2 + br~ 4 + cr~ 6 , both in three dimensions and in two dimensions, where the 
parameters a, b, and c in the potential satisfy some constraints. 
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1. Introduction 



The exact solutions to the fundamental dynamical equations play crucial roles in 
physics. It is well-known that the exact solutions to the Schròdinger equation have 
been obtained only for a few potentials, and some approximate methods are frequently 
applied to arrive at the approximate solutions. In recent years, the higher order an- 
harmonic potentials have drawn more attentions of physicists and mathematicians in 
order to partly understand a few newly discovered phenomena, such as the structural 
phase transitions [1], the polaron formation in sòlids [2], and the concept of false vacuo 
in field theory [3]. Interest in these anharmonic oscillator-like interactions comes from 
the fact that the study of the relevant Schròdinger equation, for example, in the atòmic 
and molecular physics, provides us with insight into the physical problem in question. 

For the Schròdinger equation (h = 2m = 1 for convenience) 

—V 2 ip + V(r)ip = Eip, (1) 

with the potential 

V(r) = ar 2 + br~ 4 + cr~ 6 , a > 0, c> 0, (2) 

let 

^(r,6, V )=r- 1 R e (r)Y em (0,v), (3) 

where i and E denote the angular momentum and the energy, respectively, and the 
radial wave function Ri(r) satisfies 



d 2 RÀr) 



R,(r) = 0. (4) 



dr 2 

Znojil [4,5] converted Eq. (4) into a difference equation in terms of a Laurent-series 
ansatz for the radial function 

N 

R e (r) = N r K exp[-(V^r 2 + V~cr- 2 )/2] £ h m r 2m . (5) 

m=-M 

He defined the continued fraction solutions to accelerate the convergence of the series, 
and obtained the solutions for the ground state and the first excited state. 



Kaushal and Parashar highly simplified the ansatz for calculating those solutions 

R (r) = N r K ° exp[-(v^r 2 + V~cr' 2 )/2], k = (b + 3yfi)/(2y/c), (6) 

for the ground state [6] , and 

R (r) = N ir K1 (l + (3r 2 + 7 r" 2 ) exp[-(v^r 2 + V~cr~ 2 )/2], (7) 

for the first excited state [7]. By this ansatz, the parameters in the potential (2) have 
to satisfy two constraints: 

(2v/c + ò) 2 = c[(2£ + l) 2 + 8v^c] , (8) 

and 

m [{ m - 4) 2 - 4(2«! - l) 2 ] = 64^(77, - 4), 



_ (9) 
«i = (b + 7y/c)/(2y/c), m = £(£ + 1) + 2^~c - k\ + Kl . 

where there was a sign misprint in [7] (see Eq. (13) in [7]). They set the vàlues of the 

parameters by 

£ = 0, a=1.0, c = 0.18, 6 = 0.04082, (10) 

and found that (3 = —0.1787 and 7 = 0.8485, and the energies for the ground state 
and the first excited state were E = 4.096214 and E 1 = 12.09621, respectively. Un- 
fortunately, their parameters given in Eq. (10) do not satisfy the second constraint 
(9), such that the so-called solution of the first excited state in [7] does not satisfy Eq. 
(4). As a matter of fact, they assumed that the angular momentum i is same for both 
the ground state and the first excited state, and that the normalized factor Ni 7^ 0, so 
that they must obtain, as shown in Sec. 2 of the present letter, infinite solutions for (5 
and 7 if the parameters in the potential satisfy the constraints (8) and (9). 

In our viewpoint, Kaushal and Parashar presented a good idea for studying the 
Schròdinger equation (1) with the higher order anharmonic potential (2), but their 
calculation was wrong. In the present letter, we recalculate the solutions following 
their idea, and then, generalize this method to the two-dimensional Schròdinger equa- 
tion because of the wide interest in lower-dimensional field theories recently. Besides, 
with the advent of growth technique for the realization of the semiconductor quantum 



wells, the quantum mechanics of low-dimensional systems has become a major research 
field. Almost all of the computational techniques developed for the three-dimensional 
problems have already been extended to two dimensions. 

This letter is organized as follows. In Sec. 2, we recalculate the ground state and 
the first excited state of the Schròdinger equation with this potential using an ansatz 
for the eigenfunctions. This method is applied to two dimensions in Sec. 3. The figures 
for the unnormalized radial functions of the solutions are plotted in the due sections. 

2. Ansatz 

Assume that the radial function in Eq. (3) is 



where (3 = 7 = and k — kq for the ground state, and (3 7^ 0, 7 7^ and k — K\ for 
the first excited state. Substituting Eq. (11) into Eq. (4), we have 



On the other hand, the derivative of the radial function can be calculated directly from 




(11) 



d 2 R e (r) 
dr 2 



{r 4 a(3 + r 2 [aa - (3E\ + [-aE + (3Í{Í + 1) + 7a] 
+ r~ 2 [a£(£ + l)+(3b- 7 £] + r~ 4 [ab + (3c + -y£(£ + 1)] 
+r~ 6 [ac + 7 6] +r" 8 7c}r K exp[-(v/ar 2 + y/cr~ 2 )/2]. 



(12a) 



Eq. (11), 

d 2 R e (r) 
dr 2 



{r A a(3 + r 2 [aa - (3^{2k + 5)] 
+ [-ay/a(2K + 1) + (3(2 + 3/t + k 2 - 2y/ac) + 7 a] 
+ r- 2 [a(K 2 -K- 2^/ac) + (3yfc(2K + 1) - 7v /a(2/t - 3)] 
+ r~\a^c(2K - 3) + (3c + 7(6 - 5/c - 2^fd~c + /t 2 )] 
+ r~ 6 [ac + ^^/c(2k - 7) + r~ 8 7c} r K exp[-( v /ar 2 + y/cr~ 2 )/2]. 



(126) 



Comparing the coefficients in the same power of r, we obtain 



(3[E - ^l(2k + 5)] = 0, 



(13a) 



a[E - v^(2/t + 1)] = (3[í{i + 1) + 2^[a~c - k 2 - 3/t - 2], 



(136) 



4 



a[£(£ + 1) + ly/Vc -k 2 + k}= p[-b + y/c(2n + 1)] + >y[E - y/a(2K - 3)], (13c) 

a[b - Vc(2k - 3)] = -i[í{í + 1) + - k 2 + 5k - 6], (13d) 

7[6 - Vc(2k - 7)] = 0. (13e) 
For the ground state, /3 = 7 = and « 7^ 0, we obtain a constraint (8) and 

/t = (3^ + 6)/ (2 v^), E = yoA? (6 + Vc). (14) 

For the first excited state, (3 7^ and 7 7^ 0. From Eqs. (13a) and (13e) we have 

«1 = (7y/c + b)/(2y/c), Ex = ^fc (b + 12y/c). (15) 

It is easy to see from Eqs. (8) and (15) that the right hand side of Eq. (13d) becomes 
zero, namely, a — 0. Since Kaushal and Parashar [7] assumed a^O, they must obtain 
the infmite 7 if the parameters in the potential satisfy two constraints (8) and (9). 
Now, we obtain from Eq. (13) 

ct — 0, 7 = —\Jc/a (3, (16) 

and another constraint 

b = -6<y/c. (17) 

It is easy to check that the constraints (8) and (17) coincide with the constraints (8) 
and (9). 

Setting í = and a = 1.0 for comparison with Znojil [4] and Kaushal-Parashar [7], 
we obtain 

b = -11.25, Jc = 1.875, 7 = -1.875/5, 

v 1 (18) 

/t = -1.5, /ti = 0.5, E = -2, E 1 = 6. 

Thus, the radial functions -Ro^(r) for the ground state and fí^\r) for the first excited 
state are 

R^\r) = iVor- L5 exp{-(r 2 + 1.875r" 2 )/2}, 

i?W(r) = Ar ir - - 5 (r 2 - 1.875r- 2 )exp{-(r 2 + 1.875r" 2 )/2}, 



where the normalized factors are calculated by the normalized condition: 

í°° \R^(r)\ 2 dr = 1, i = and 1. (20) 
Jo 

Without loss of any main property, we show the unnormalized radial functions in Fig. 
1 and Fig. 2. 

Furthermore, if the angular momentum £' for the first excited state is different from 
the angular momentum £ for the ground state, equation (16) and the constraint (17) 
become 



P = Aay/E/[e'(e' + 1) - 1(1 + 1) - 4(6 + 6y/c)/y/^, 

1 = Aa y /c/[l'{e + l)-£(£ + l)}, 

[£'(£' + 1) - 1(1 + 1) - 2(6 + 4^)/^ /(32^/Ec) 

= [£'(£' + 1) - £(£ + 1) - 4(6 + Gv^/v^ 1 + [f (f + 1) - 1(1 + l)]" 1 . 

Setting a = 1.0, £ = and £' = 1, we obtain 

6 = -4.2011, c = 0.75878, k = -0.91144, Kl = 1.08856, 
P = -1.47683 a, 7 = 1.74216 a, E = -0.82288, £1 = 7.17713. 



(21) 



(22) 



3. Solutions in two dimensions 
For the Schròdinger equation in two dimensions with the potential, 

V(p) = ap 2 + bp' 4 + ep" 6 , a > 0, c>0, (23) 

let 

^(p, <p) = p- 1 / 2 R m (p)e ±i ^, m = 0, 1, 2, • • • , (24) 
where the radial function R m (p) satisfies the radial equation 



d 2 R m (p) 



+ 



E - V(r) - 



m 2 - 1/4 



i?m(p) = 0. (25) 



dp 2 

Making the ansatz for the radial functions of the ground state and the first excited 
state: 

i4°)(p) = N o/ f° expHvV + v^p~ 2 )/2], 

(26) 

Rg\p) = N lP l* (a + p 2 + 7P" 2 ) expHvV + v^P" 2 )/2], 



where 7 7^ 0, and substituting Eq. (26) into Eq. (25), we have 
P[E - y/a(2 K + 5)] = 0, 

a[E - y/a{2K + 1)] = f3[m 2 - 1/4 + 2y/iïc~ - k 2 - 3/t - 2], 

a[m 2 - 1/4 + 2^ - k 2 + «] = /?[-& + ^(2/t + 1)] + 7 [£ - ^(2/t - 3)], (27) 
a[& - v/c(2/t - 3)] = -7[m 2 - 1/4 + 2^/ac - k 2 + 5/t - 6], 
7 [6 - v^(2« - 7)] = 0. 

Hence, if the angular momentum m of the ground state is the same as that of the 
first excited state, we obtain from Eq. (27) 



{2^l· + bf = Ac[m 2 + 2^\, b = -6y/c, 
K = (3yft + b)/(2y/c), E = ^fc (b + Ay/c) , 

ki = (7Vc + 6)/(2Vc), E 1 = ^~[c (b + 12y/c), 

a = 0, 7 = -y/c, 

If m — and a = 1.0, the vàlues of the corresponding parameters are 

b = -12, c = 4, 7 = -2, 

Kq = —1.5, «i = 0.5, Eq = —2, Ei = 6, 



(28) 



(29) 



The unnormalized radial functions are shown in Fig. 3 and Fig. 4. 

To summarize, we discuss the ground state and the first excited state for the 
Schròdinger equation with the potential V(r) = ar 2 + br~ A + cr" 6 using a simple 
ansatz for the eigenfunctions. Two constraints on the parameters in the potential are 
arrived at from the compared equations. This simple and intuitive method can be 
generalized to the other potentials, such as the sextic potential, the octic potential, 
and the inverse potential. 
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